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ABSTRACT 

Fomalhaut and e Eridani are two young, nearby stars that possess extended 
debris disks whose structures suggest the presence of perturbing planetary ob- 
jects. With its high sensitivity and stable point spread function, Spitzer/IRAC 
is uniquely capable of detecting cool, Jupiter-like planetary companions whose 
peak emission is predicted to occur near 4.5 /xm. We report on deep IRAC imag- 
ing of these two stars, taken at 3.6 and 4.5 fim using subarray mode and in all 
four channels in wider-field full array mode. Observations acquired at two differ- 
ent telescope roll angles allowed faint surrounding objects to be separated from 
the stellar diffraction pattern. No companion candidates were detected at the 
reported position of Fomalhaut b with 3a model-dependent mass upper limits 
of 3 Mj (for an age of 200 Myr). Around e Eridani we instead set a limit of 4 
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and < 1 Mj (1 Gyr model age) at the inner and outer edge of the sub-milhmeter 
debris ring, respectively. These results are consistent with non-detections in re- 
cent near-infrared imaging searches, and set the strongest limits to date on the 
presence of planets outside e Eridani sub-millimeter ring. 

Subject headings: infrared: stars - planetery systems - stars: individual (Fomal- 
haut, e Eridani) 



Introduction 



Debris disks are dust clouds in extrasolar planetary systems produced by ongoing col- 
lisions between small bodies analogous to asteroids and comets. Their internal structure is 
of considerable scientific interest, for it can directly betray the presence of planets. Narrow 
rings, warps, central holes, and azimuthal asymme tries can be related to the likely mas s 
and orbital characteristics of the perturbing object (ILiou &: Zooklll999l : IWyatt et al.lll999l ). 
Nearby debris disks are thus ideal targets for deep imaging searches for extrasolar planets. 

Fomalhaut {a PsA; IRAS 22549- 2953) is a bright A3 V star 7.7 pc distant, with an esti- 
mated age between 100 and 300 Myr (iBarrado y Navascueslll998l ). An inclined circumstellar 
dust ring, with an inner and out er edge diameter of 260 and 310 AU respectively, was first 
resolved bv [Holland et al.l (|l998f) at 850 yum. A modest ring brightness asymmetry seen at 
450 /im ( Holland et al. 20031 ) be comes increasingly prom inent toward shorter wavelengths in 
Spitzer Space Telescope images (IStapelfeldt et al.ll2004j ). consistent with models of an eccen- 
tric ring that is warmed at p eriastron. This ni odel was subsequently confirmed in detail by 
HST scattered light imaging (IKalas et al.ll2005l ). The sharpness of th e ring iiiner ed ge is con- 
sistent with dynamical sculpti ng by a giant planet orbiting nearby (jQuillenll2006l ). recently 
detected by iKalas et al.l (120081 ) at optical wavelengths. 



e Eridai ii (IRAS 03305-0937 ) is a K2V star only 3.2 pc distant, with an estimated age of 
~ 850 Myr (jPi Folco et al.ll2004l ). A su bmillimeter dust ring , 220 AU in diameter and pro- 
jected close to face-on, was discovered by lGreaves et al.l (119981 ). Azimuthal clumps seen along 
the ring have been suggested to be dust trapped in orbital resonances with a giant planet 
( lOzernoy et al.l l2000l : iQuillen &: Thorndikd 120021). but the brightness and location of these 



clumps have not been reproducible (jSchiitz et al.ll2004l : iGreaves et al.ll2005l : iBackman et al. 



20091 ). Far- infrared imaging and spectroscopy with Spitzer indicate the presence of two inte- 
rior dust rings, with additio nal planets suggested to clear the radial gaps in the dust distribu- 



searches ( iMacintosh et al. 



tion (iBackman et al.l 12009). The star has been the target of many near-in f rared companion 



20031 : ICarson et all l2005l : iNakaiima et al.ll2005l : lltoh et all 12006 



- 3 - 



Biller et al.ll2007l : Ijanson et al.l 120071 . 12OO8I : iHeinze et al.ll2008l : iKenworthv et al.ll2009f ). Sev- 
eral of these stu dies focused on det ecting the radial velocity planet that that orbits about 
1" from the star (IHatzes et al.ll2000l ). but without success. 



Models for the emergent spectra of young giant planets show a prominent ernission 
peak near 4.5 Atm, caused by a gap in CH4, H2O and NH3 opacity (IBurrows et al.l |2004J : 



Fortney et al.ll2008l ). Spitzer /IR AC is well-suited to searching for such objects, which in 
broad-band photom etry will appear bright in channel 2 and much fainter in channel 1 



(IPatten et al 



2006). T dwarf companions to two nearby stars were discovered in this way 



by iLuhman et al.l (2007). A deep IR AQ search for planetary companions to e Eridani was 
conducted by iMarengo et al.l (120061 ). No companion candidates were found to a model- 
dependent mass limit of 1 Jupiter mass, in the region r > 25" (80 AU). Inside this radius, 
the imaging sensitivity was compromised by bright star artifacts. 

The IRAC subarray imaging mode offers rapid readouts and exposure times as short 
as 0.02 s. To probe the 5-20" region where the Fomalhaut and e Eridani debris rings are 
found (and where planetary companions perturbing the rings are expected), we conducted 
a program of deep IRAC subarray imaging of these two stars. In this paper, we present 
the results of our IRAC subarray imaging search, a new analysis of the wider-field IRAC 
full array dataset, limits to the surface brightness of the debris disks at IRAC wavelengths, 
and limits to the mid-infrared brightness of Fomalhaut b and of planets in proximity of the 
e Eridani rings. 



2. Observations and data reduction 



The observations are summarized in Table [H The IRAC subarray observations were 
carried out as a part of Spitzer General Observer (GO) program 30754 using only band 1 
(3.6 fim) and band 2 (4.5 fim). Each observation requested 255 repeated stacks of 64 x 0.02 s 
frames (0.01 s integration time) on 9 dither positions, spaced over the subarray on a small- 
scale Reuleaux triangle. The choice of the frame time and dither pattern was guided by the 
need to reduce saturation as much as possible and provide the best image sampling on the 
under-sampled IRAC pixel grid, in order to allow the removal of the stellar Point Spread 
Function (PSF). Fomalhaut was observed with this sequence on 2006 November 24 and 2006 
December 26, with a 11.67° clockwise roll angle offset between the two epochs, e Eridani was 
observed on 2007 February 16 and 2007 September 13, providing a 169.03° roll angle offset. 
The total on-source integration time for each sequence of 9 x 225 x 64 = 146, 880 frames was 
1468.8 s. The first Fomalhaut sequence executed on 2006 November 24, however, was cut 
short due to a bug in the Spitzer time allocation software, skipping the last 77 repeats in 
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the final dither position while observing in band 2. The integration time for the Fomalhaut 
first epoch 4.5 fim observation is thus only 1,419.5 s. The total field of view for the subarray 
observations is ~ 44". 

IRAC full array observations of Fomalhaut took place on 2004 Oct 31 and 2004 Nov 28, 
as a part of Guaranteed Time Observer (GTO) program 90. The roll angle offset between 
the two epochs was 11.03° clockwise. At each epoch, the observations consisted of four 12 s 
frames (10.4 s integration time) taken at each of 12 dither positions spaced over the array 
on a small-scale Reuleaux triangle. The total integration time was 499.2 s in each of the 
four IRAC bands. The total field of view imaged in each band was ~ 5.78'. e Eridani was 
observed on 2004 Jan 9 and 2004 Feb 17 with a roll angle offset of 20.15° clockwise. The 
total integration was of 3369.6 sec in each of the 4 IRAC bands, obtained using individual 
12 sec frame-time observations on a 36 points Reuleaux small scale dither pattern, repeated 
9 times in each position . The e Eridani full frame IRAC data has been previously presented 



m 



Marengo et al.l (120061 ). We include here a new analysis of this data to take advantage of 
the latest photometric calibration of IRAC, and to ensure a consistent data reduction of the 
full frame images with the subarray observations. 

The two stars, in both subarray and full frame observations, are saturated. While the 
full frame data are saturated within a radius of ~ 6" from the star, the much shorter frame 
time of the subarray observations restricts saturation to the central 1-2 IRAC pixels (~ 1.5- 
3" radius). This allows to search for low mass companions in the subarray data at closer 
distances from the star than in the full frame data. 



2.1. Data Reduction 

Basic data reduction for all the observations was performed with the Spitzer Space 
Center (SSC) IRAC Pipeline version S15, which produced Basic Calibrated Data (BCD) 
frames and data quality masks for each individual full fr ame and subarray ex posure. The 



full frame data were reduced using the IRACproc package (ISchuster et al.ll2006l ). to obtain a 
single flux calibrated mosaic combining all the individual exposures in each epoch and each 
band, on a pixel grid with 0.24"/pixel resolution. IRACproc is based on the SSC mosaic 
software MOPEX and provides enhanced outlier (cosmic rays) rejection. 

We reduced the subarray data with our custom software, due to the very large number of 
frames acquired for each exposure, and to mitigate strong "muxstripe" artifacts in the BCDs. 
The muxstripe artifacts appear as a "jailbar" pattern repeating with a 4 columns cadence, 
and is triggered in the InSb arrays (bands 1 and 2) by the flux of bright sources unbalancing 
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the 4 multiplexer readout channels. The striping pattern spans the whole array, but is 
stronger below the triggering source, where an exponential decay pattern is superimposed 
with the same 4 column cadence. The muxstripe is poorly characterized (does not scale 
linearly with the source brightness) and its strength depends on the Fowler sampling numbers 
of the exposures. We unexpectedly found that for the brightness of our sources and the chosen 
frame time the muxstripe artifacts were dominating all other sources of noise (as much as 
~ 50 times above the expected read-out noise in the final mosaics). 

An example of the problem is shown in Figure [1] for a stack of 255 x 64 co-added (with 
a basic temporal outlier filter) frames (panel a). The first step we took to characterize 
the artifacts was to isolate them by subtracting the stellar PSF, which we constructed by 
combining the subarray observations of 25 bright Cepheid stars obtained as part of the 
General Observer program 30666. These Cepheids are bright enough to provide a high 
S/N PSF without being saturated. Observed with frame times of 0.1 and 0.4 s, they still 
suffer from muxstriping, but at a much lower level than our 0.02 s frames, and without the 
presence of the exponential decay. The "Cepheid PSF" was aligned with each 255 x 64 co- 
added image using the diffraction spikes (typically with a precision of better than ~ 0.2 ") 
and subtracted. Figured] (panel b) shows the same frame after PSF subtraction. While the 
PSF structures are not fully removed (due to the different observing parameters and source 
colors), the muxstriping is more clearly seen. To characterize this artifact we first equalized 
the column offset pedestal between the 4 multiplexer readout channels. We then fitted each 
column separately, in the lower part of the image, with either an exponential function of 
the form y = b ■ exp(— s/a) -|- c or a linear slope y = bx + c {in case of not converging 
exponential fit). The result of this operation is a "correction matrix" (panel c in Fig[T]) that 
is then applied to the PSF-subtracted 255 x 64 co-added stack (panel d). The correction 
was derived independently for each of the 9 dither positions, as the artifacts depend on the 
position of the star on the pixel grid. 

The 9 individual images obtained for each dither position with the procedure described 
above were then aligned (using the centroids derived when removing the "Cepheid PSF") 
and co-added to produce one image per band, and per epoch, for each source. 

2.2. PSF roll-subtraction 

We have combined the data for the two epochs in one single image for the full frame 
and subarray observations, separately. By taking advantage of the different roll angle, we 
have removed the stellar PSF while preserving the other point sources in the field. 
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In the case of the subarray data we ahgned the two images obtained for each epoch. 
Even though these images have been aheady PSF-subtracted with the "Cepheid PSF" , this 
subtraction leaves significant residuals because of the different dither pattern adopted. To 
remove these residuals, we subtracted the aligned images one with the other. After masking 
in each image the row centered with the subtracted star (where the muxstripe correction is 
less accurate), we then rotated the subtracted images to have N up and E left. We then 
averaged them to obtain one single image for each source, in each band (Figure [2]). Because 
of this procedure, point sources in the field of these images would appear as one bright spot 
sided by two negative sources having half height: inspection of the images in Figure [2] clearly 
shows that we do not detect any source within the subarray field of view towards either stars. 

In the case of the full frame images, we have taken advantage of the availability of 
a high S/N IRAC PSF to improve the PSF removal beyond what is achievable with the 
standard 2 epochs roll-subtraction technique. We derived the full frame IRAC PSF by 
combining observations of all the stars in the GTO program 90 (2 epoch images of e Eridani, 
Fomalhaut and Vega, 1 epoch image of e Indi) with the exclusion of (3 Pic (because of the 
presence of its debris disk, detected at IRAC wavelengths), plus an image of Sirius obtained 
as part of the IRAC calibration program. We first reduced the individual PSF star images 
using IRACproc, producing mosaics on a pixel grid rotated as the IRAC arrays (to preserve 
the orientation of the PSF features) with a pixel scale of ~ 0.24"/pixel. We masked any pixel 
with flux higher than 80% of the IRAC saturation limit, to ensure linearity. We then aligned 
all the images together, using the diffraction spikes as reference (this typically provides an 
accuracy of ~ 1/10 of the pixel grid, or ~ 0.02"). Finally, we rescaled all the images to the 
same reference (one of the images of Vega) and then co-added them together using a median 
filter to reject all point sources in the field. This produced a single PSF image in each band, 
calibrated to have the same surface brightness (in MJy/sr) as Vega, cleaned of all point 
sources in the field. Note that, as a byproduct of this procedure, the normalization factor 
of the PSF stars are nothing else than the flux ratios of these stars with Vega. These ratios 
can be express as a measure of the Vega magnitudes of the stars, completely independent 
from the absolute IRAC calibration (because they are the ratios of IRAC images of the stars 
with an actual image of Vega), and are listed in Table [H The typical accuracy of these 
magnitudes is ~ 0.01 mag. These PSFs are available at the SSC web sit^. 



The procedure described above is similar to the one adopted in iMarengo et al.l (120061 ). 
but results in a more accurate PSF because of the addition of Sirius and because of a better 
masking procedure, preserving linearity closer to the center (ensuring more accurate PSF 



^ http :/ /ssc . spitzer .caltech.edu/ir ac/psf. html 
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subtractions near the star, which is crucial for this program) . While this PSF has a very high 
S /N and is cleaned of the point sources in the field, it still suffer from uncorrected electronic 
artifacts (column pull-down, muxbleed, bandwidth effect and banding). These artifacts are 
not linearly dependent on the source flux, and are not removed by simple PSF-subtraction. 
To improve this situation, we have derived two specialized version of the PSF, suitable to 
remove both electronic and optical PSF features, one optimized for Fomalhaut, and the other 
for e Eridani. 

To achieve this, we first aligned and rescaled the two images of each source (one for each 
epoch, left in the IRAC array orientation, so that the structures of the PSF will overlap), 
with the PSF. We then median combined the three-image stack, obtaining a new PSF that: 
(1) in areas free of electronic artifacts, it is the median of the original PSF with the two 
epoch images; and (2) where electronic artifacts occur, it samples values of the artifacts from 
one of the two source images. This optimized PSF is free of background sources (which are 
median filtered) and at the same time has the electronic artifacts scaled appropriately for 
the source fluency. We then aligned and subtracted the optimized PSF from each of the 
two epochs, achieving a subtraction of both the optical and electronic structures of the PSF. 
This subtraction is of higher quality that the standard two epochs roll-subtraction, because 
the optimized PSFs are free from background sources. 

After subtracting each of the two epoch images from their optimized PSF, we projected 
the result to have N up and E left, and then averaged the rotated images. A sample PSF 
roll-subtracted image is shown in Figure [3] (bottom panel, the top panel shows the image 
before PSF subtraction) for Fomalhaut at 4.5 /xm. Figures of the final PSF roll-subtracted 
full frame images for both sources in all bands are available in the electronic version of this 
paper (Plates 1 to 8). 



3. Results 



Figure [2] and [3] sho w that we do not detect any point source at the location of Fo- 
malhaut b indicated by iKalas et al.l (120081 ). We can however use our PSF roll-subtracted 
subarray and full frame images to set limits on the brightness of the planet in the IRAC 
bands, and of point sources and extended emission around Fomalhaut and e Eridani. 
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3.1. Sensitivity 

Figure m and [5] show the 3a radial sensitivity for our subarray and full frame images. We 
estimated the sensitivity curves by measuring the root mean square (RMS) noise in circular 
annuli having a width dr ~ 2FWHM(A), where FWHM(A) is the Full Width Half Maximum 
radius of the IRAC PSF in each band. From the RMS noise (which is in surface brightness 
units of MJy/sr) we have then derived the limiting magnitude in each band, and for each 
annulus, as follows: 



mi 



im 



-2.5 log 



Svr (FWHM/2)^ • RMS 
Fo(r = FWHM/2) 



where Fq is the flux of the IRAC PSF, normalized as Vega, within a circular aperture 
with the diameter of the PSF FWHM, and the factor 3 has been introduced to obtain 
3(7 sensitivities. For a Gaussian noise pattern, these 3a sensitivities would guarantee a 
99.7% detection probability. Given that the PSF-subtraction residual noise is not Gaussian, 
however, we have tested these sensitivity curves by "planting" point sources of different 
brightness in the PSF roll-subtracted images, verifying that the plotted 3a curves indeed 
guarantee the detection of the sources in all cases. 

The two figures also show the expected magnitudes o f 1, 2, 5 and 10 Mj planets of 
200 Myr and 1 Gyr of age, interpolated from models by iBurrows et al. Similar 
estimates can be obtained using the more recent models by iFortney et al.l (120081). The 



vertical line in Figure H] indicates the projected separation of Fomalhaut b from lKalas et al. 



([20081), while the vertical lines in Figure [5] sho w the inner and outer radii of the e Eridani 
submillimeter ring (from iBackman et al.l 120091 ). 



The plots show that the subarray images are sensitive to point sources as close as ~ 3" 
from the central star. The sensitivity at such small radii is generally poor (~ 10 mag for 
Fomalhaut and ~ 11 mag for e Eridani), not sufficient to detect planetary mass bodies around 
either stars, but enough to detect T dwarfs. As explained in Section [2] the subarray images 
sensitivity is limited by the residual errors from the muxstripe correction, and "fiattens-out" 
at radii larger than 15". The extent of the muxstripe artifacts in subarray mode was not 
expected to be so strong at the time these observations were proposed. For radii larger than 
~ 6", however, our new analysis of the full frame PSF roll-subtracted images provides a 
better sensitivity. Between ~ 6 to 50" from the star, our full frame images are limited by 
the PSF subtraction noise, while for radii larger than ~ 50" the sensitivity is limited by the 
noise in the background. 
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The 3(7 sensitivity limits at the position of Fomalhaut b (Figure H]) show that the planet 
has a brightness lower than 0. 5, 0.7, 2.05 and 1.39 m Jy at 3.6, 4.5, 5.8 and 8.0 fim respectively 
(see Table S}. According to iFortney et al.l (120081 ) this implies a model-depend ent mass of 
< 3 M j (for a ~ 200 Myr age), in agreement with a similar limit inferred by iKalas et al. 



( I2OO8I ) from ground-based near-IR data using the same models. Our images provide the first 
brightness upper limits for the planet at 4.5, 5.8, and 8.0 /zm (an L' band limit is given in 
Kalas et al.ll2008h . 



From the noise RMS derived above, we can also set limits to the surface brightness of the 
debris disks around both stars, listed in Table [51 Note how our sensitivity at the location 
of the Fomalhaut ring NW ansa is as much as 5 times better (1.7 magnitudes) than the 
sensitivity at the location of the SE ansa. This is because the SE ansa happens to overlap 
with high noise residuals from PSF "pull-down" electronic artifacts, locally decreasing our 
sensitivity below the circular average level plotted in Figure HI Table [5] also shows our 
sensitivity limits for the e Eridani second asteroid belt (at 20 AU radius) and for the inner 
and outer edge of th e sub-millimeter ring (at 35 and 100 AU respectively) described by 
Backman et^ fl2009l ). 



Macintosh et al.l (120031 ) and iJanson et al.l (120081 ) derive limits for the mass of planetary 
mass bodies inside the inner rim of the e Eridani sub-millimeter ring of 5 and 3 Mj re- 
spectively, compared to our limit M > 4 Mj (mass limits estimated for ~ 1 Gyr models). 
Outside the sub-millimeter ring our limit (M > 1 Mj, 1 Gyr model age) is s u perior to any 



other available observation, including our previous analysis in iMarengo et al 



3.2. Point Source Photometry 

To search for low mass companions around Fomalhaut and e Eridani we have measured 
the photometry of all point sources detected in the PSF roll-subtracted images of both stars. 
As shown in Figure H] and E] we are sensitive to substellar objects with mass lower than 
1 Mj outside the debris rings of both stars (from models with age of 200 Myr and 1 Gyr 
respectively). 

To measure the photometry of all the point sources in the field of view of the full 
frame images of both stars, we first converted the data from units of surface brightness 



^Note that Figure 7 in lMarengo et al.l ( 20061 ) shows the sensitivity in an area of the PSF-subtracted image 
far from diffraction spikes or other artifacts. Figure[5]in the present work shows instead a circularly averaged 
sensitivity, and is superior to the equivalent measure performed on our 2006 PSF-subtracted images. 
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flux density to magnitudes, by dividing by the flux conversion factors found in tlie image 
fieaders of tlie BCD, and by multiplying by the effective integration time. The IRAC absolute 
calibration is based on observations of standard stars measured with aperture photometry 
using a source aperture with a radius of 10 native IRAC pixels (~ 12.2") in each channel 
(this is equivalent to state that the aperture correction of 10 IRAC pixel radius apertures 
is exactly 1 according to this convention). The background was estimated using an annulus 
centered on the so urce position with an inner radius and width of 10 IRAC native pixels 
(IReach et al.ll2005l ). Because our two targets are located in semi-crowded fields, we chose 
to use a smaller source aperture with a radius of 3.6" in order to avoid contaminating flux 
from other, nearby sources. An additional benefit of using a smaller source aperture is an 
improved S/N for many of the fainter sources. For the background estimation, we used 
the annulus with an inner radius of 4.8" and a width of 2.4". Photometry was extracted 
using the IRAF phot package. We determined an aperture correction for these nonstandard 
aperture sizes by comparing the photometry from a set of bright stars in our images using the 
nominal parameters. We then determi ned the neces s ary z ero points to produce photometry 
in the standard calibration defined by iReach et al.l (120051 ). The zero points used are listed 
in Table O 

We have first analyzed the colors of all sources detected in multiple IRAC bands. Given 
the lower sensitivity at 5.8 and 8.0 /xm, and the CH4 absorption expected at 3.6 fj,m for T- 
dwarfs and gas giant planets, this limits our discovery space to objects with a 200 Myr model 
dependent mass of ~ 2 Mj around Fomalhaut, and ~ 4 Mj around e Eridani (1 Gyr model 
age). We used the same /c-Nearest Neighbor method search technique base d on the colors 



and a bsolute magnitudes of b rown dwarf and planet mo dels we adopted in iMarengo et al. 



(l2006l ). described in detail in iMarengo fc Sanchea (120091 ). None of the sources detected in 
all four bands of full array imaging had the colors and absolute magnitudes expected for 
sub-stellar mass objects. 

All sources that were de tected at 3.6 and 4.5 u m, b ut missed at 5 . 8 and / or 8.0 /im, have 
also been rejected. Based on lBurrows et al.l (120031 ) and lFortney et al.l (120081 ) models none of 
them possess a [3.6] — [4.5] color red enough to be a planetary mass object for their measured 
4.5 fim magnitude. 

We then focused on sources detected at 4.5 / im but not detect e d at 3.6 Atm: these 
sources can be either very red extragalactic sources (IHuang et al.ll2004j : iBarmby et al.ll2004l ) 
or methane dwarfs and planets. Red, mass losing bac kground giants (such as Asymptotic 
Giant Branch stars) have positive [4.5] — [8.0] colors (IMarengo et al.l l2008l ). and are thus 
distinguishable from T dwarfs and planets that have instead a 8.0 /xm flux equal or lower 
than the 4.5 micron flux. Table [6] and [7] lists all 3.6 /xm "dropout" sources for Fomalhaut 
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and e Eridani respectively. Note how several of them can be discarded because their flux 
at 5.8 and 8.0 fim makes them likely to be background mass losing giants. Comparison 
with other datasets could enable the rejection of other dropout sources based on (lack of) 
common proper motion with the primary star. To our knowledge, however, none of the 
dropout sources in Table [6] and [7] has been detected in any other observation at optical or 
infrared wavelength. 

Figure E] shows the [3.6] — [4.5] vs. [4.5] color-magnitude diagram of all sources detected 
in both IRAC bands 1 and 2 (data points with error-bars) and of all 3.6 /im dropout sources 
(arrows). The [3.6] — [4.5] color of the dropout sources is a lower limit, estimated using the 
local 3.6 /xm 3a sensitivity derived in Section [3TT1 The dropout sources within the Fomalhaut 
field of view h a .ve 4. 5 /xm magnitudes expected for < 1 Mj 200 Gyr planets, according to 
Burrows et al.l (120031 ) models. However, th eir [3.6] — [4.5] col ors are still compatible with the 
colors of red extragalactic sources (e.g. see lStern et al.ll2007l ). Similarly, the dropout sources 
in the e Eridani field (that, if they were ~ 1 Gyr planets, would have a mass < 2 Mj) also 
have color limits compatible with being background galaxies. 

Figure [7] shows the position of the 3.6 fim dropout sources within 50x50" from Fo- 
malhaut and e Eridani, and their relative position with respect to the debris disks. This 
area corresponds to ~ 385x385 AU at the distance of Fomalhaut, and ~ 160x160 AU at 
the distance of e Eridani. Three dropout sources are detected within 300-400 AU from the 
Fomalhaut ring, but none of them is in close proximity, or inside, the debris ring. Of these 
three sources, one has [4.5] — [5.8] and [4.5] — [8.0] color larger than ~ 1.5, which is charac- 
teristic of background mass losing giants. The oth er two sources are c ompanion candidates 
with mass lower than 1 Mj (according to 200 Myr lBurrows et al.ll2003l models). Of the four 
dropout sources detected near e Eridani, one is inside the sub-millimeter ring. This source, 
as well as two of the sources detected within ~ 100 AU from the outer radius of the e Eridani 
sub-millimeter ring, is likely a mass losing giant in the background. The remaining source, 
~ 90 AU from the out er rim of the debris disk, is a planet candidate with a mass as low as 
2 Mj (based on 1 Gyr lBurrows et al.l l2003l models). 



Of all the dropout sources listed in TableO only 5 were detected at 4.5 fim in lMarengo et al. 



( 120061 ) ■ due to the lower quality of the PSF su btraction in that work. All the other 3.6 /im 
dropout sources found in lMarengo et al.l (120061 ) . Table 3, have been detected above 3a in our 
current PSF roll-subtracted images, and have been excluded as planet can didates. Thanks 



to our better sensitivity at 4.5 fim, we detect 3 sources that were missed by lMacintosh et al. 



( 120031 ) in their field of view. Of these sources, two are 3.6 /im dropouts (shown in Figure [TJ E 
of e Eridani). The colors of these sources suggest that they are background stars or galaxies. 
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Discussion and Conclusions 



The phys ical origin of t h e opt ical light detected by iKalas et al.l (120081 ) from Fomalhaut 
b is unclear. iFortney et al.l (120081 ) model of the thermal emission from a 400 K object of 
~ 1.2 Jupiter radius (corresponding to a 3 MJ object at age 200 Myr) can reproduce the 
0.8 fim flux density observed by HST. This model is shown in Figure [HI alo ng with the 



Spitze r mid-infrared upper limits from this study and the photometry from iKalas et al. 



(120081 ). As already noted by the latter authors, this model predicts a 1.6 /im brightness 
~ 5x larger than the observed upper limit, and a 0.6 /im flux density orders of magnitude 
fainter than the observed HST detection (even taking into account the observed 0.6 /im 
variability, possibly related to variable Ha emission from a hot planetary chromosphere). 
Reconciling the thermal emission model to the observations would require an additional 
physical mechanism to produce the 0.6 /im emission, and that atmosphere models be revised 
to suppress the expected 1.6 /im emission. The Spitzer 4.5 /im upper limit lies on top 
of this thermal emission model, and thus strengthens these conclusions. Any revision of 
the thermal emission model to account for the 1.6 /im non-detection is constrained by our 
results. In particular, our 4.5 /im upper limit does not allow a lot of room for the suppressed 
1.6 /im luminosity to emerge instead through the largest low-opacity spectral window in 
methane-dominated atmospheres. The strength of this constraint can only be evaluated 
through new model atmosphere work directed toward finding a thermal emission solution to 
the properties of Fomalhaut b. In the interim, reflection from a circumplanetary dust disk 
remains the simplest model to explain the observed fluxes from the object. 

Our surface brightness upper limits (Table [5l) for the Fo n aalhau t ring can be compared 
to the contrast seen in optical scattered light by iKalas et al.l ( 120081 ). On the NW ansa, our 
3.6 /im limiting surface brightness (in mag/arcsec^) is 17.5 mag fainter than the star itself. 
At 0.6 /im, the HST detected surface brightness of 21.0 mag/arcsec^ is 19.9 mag fainter than 
the star. To escape detection, the V-L' color of the ring must therefore be less than 2.4 mag. 
V-L' colors of debris disks have however still not been measured, so the strength of our color 
constraint is unclear. However, t he reddest class of Kuip er Belt objects in our solar system 
has V-K colors less than 1.5 mag ( Cruikshank et al.ll2007l ). If the Fomalhaut dust has similar 



properties, it thus should not have been detected in our IRAC images. 

Similarly, the expecte d surface brightness of the e Eridani sub-millimeter ring, as derived 
by iBackman et al.l (120091 ) , is at least 1 or 2 o r ders of magnitude below our sensitivity, in 
agreement with the 0.011 MJv/sr lProffitt et al.l (120041) li mit se t by STIS camera observations 
in scattered light. According to the iBackman et al.l (120091 ) model, the e Eridani 20 AU 
asteroid belt is expected to have a surface brightness of ~ 21.7 mag/arcsec^ (~ 0.3 MJy/sr) 
in the V band, corresponding to ~ 0.02 MJy/sr at 3.6 /im. This is also well below the 
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sensitivity limits listed in Table [5] and cannot be detected in our IRAC images. 

None of the 3.6 /im "dropout" sources listed in Tableland [7] have colors and magnitudes 
strongly suggesting that they are planetary mass companions of Fomalhaut or e Eridani, 
rather than red background objects (even though repeated observations aimed to measure 
common proper motion are required to clarify the issue). Our detection limits for point 
sources outside the debris disks of Fomalhaut and e Eridani thus imply the likely absence of 
an y widely separated c ompanions of the two stars with mass larger than 1 Mj (according 
to iBurrows et al.l l2003l models of 200 Myr and 1 Gyr age, respectively). Inside the rings 
(and down to the 3-6" saturation radius of our images) our limits obtained with the Spitzer 
85 cm aperture telescope are comparable, or superior, to model-dependent mass detection 
limits from infrared observations obtained at 8-10 m class ground-based telescopes. 



We thank John Krist and James Graham for their advice on PSF subtractions and 
the application of the Fortney atmosphere models to our dataset. This work is based on 
observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology under a contract with NASA. Support for this 
work was provided by NASA under Spitzer General Observer grants 30754 to the Harvard- 
Smithsonian Center for Astrophysics and to JPL. 
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Table 1. Observing Log 



Source 


Mode 


AOR KEY 


Obs. Date 


Tot. Exp. 


Wavelen^ 


^ths 








[JD-2400000.5] 


[sec] 


[lira] 




Fomalhaut 


Full 


4875776 


53337.306 


499.2 


3.6, 4.5, 5.1 


3, 8.0 


Fomalhaut 


Full 


9015040 


53309.253 


499.2 


3.6, 4.5, 5.! 


3, 8.0 


Fomalhaut 


Subarray 


18951936 


54063.253 


1468.8 


3.6, 4.5 




Fomalhaut 


Subarray 


18952192 


54095.207 


1468.8 


3.6, 4.5 




e Eridani 


Full 


4876032 


53013.836 


3369.6 


3.6, 4.5, 5.i 


3, 8.0 


e Eridani 


Full 


4876288 


53052.781 


3369.6 


3.6, 4.5, 5.! 


3, 8.0 


e Eridani 


Subarray 


18951424 


54147.336 


1468.8 


3.6, 4.5 




e Eridani 


Subarray 


18951680 


54356.199 


1468.8 


3.6, 4.5 





^Total Exposure time for the 4.5 /um image is 1419.5 s 
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Table 2. Source PSF-fitting Photometry 



Source 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


Vega^ 


0.00 


0.00 


0.00 


0.00 


Fomalhaut 


0.98 


0.98 


0.98 


0.98 


e Eridani 


1.60 


1.63 


1.61 


1.59 


e Indi 


2.10 


2.15 


2.12 


2.07 


Sirius 


-1.38 


-1.38 


-1.38 


-1.35 



Used as reference, it has Vega- magnitude by definition 
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Table 3. IRAC Photometric Calibration^ 







IRAC band 




Item 


3.6 /im 


4.5 /im 


5.8 /im 


8.0 /im 


Isophotal A [/xm] 


3.550 


4.493 


5.731 


7.872 


FLUXCONV [(MJy/sr)/(DN/s)] 


0.1088 


0.1388 


0.5952 


0.2021 


GAIN [e/DN] 


3.3 


3.71 


3.8 


3.8 


Zero point magnitudes^ 


16.981 


16.512 


16.013 


15.439 




280.9 


179.7 


115.0 


61.1 



^Based on the IRAC Data Handbook ver. 3.0 (2006) 

^For pixel size 0.24"/pix, include aperture correction for 3.6" aperture and sky annulus 
with 4.8" and 7.2" inner and outer radii 
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Table 4. 3a limiting Sensitivity at the Fomalhaut b Radius 



Mode 


Fiim{3.6iim) 


Fiim{4:.5iim) 


Fiim{5.8iim) 


Fiim{S.OfJ.m) 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


Full frame 


0.50 


0.70 


2.05 


1.39 


Subarray 


2.54 


1.83 
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Table 5. Fomalhaut and e Eridani Disks brightness limits 



Mode 




Siimi4..5iJin) 


Siim{5.Snin) 


SHm{8.0iJ.m) 




[MJy/sr] 


[MJy/sr] 


[MJy/sr] 


[MJy/sr] 


Fomalhaut ring (NW ansa) 


1.19 


1.43 


3.84 


2.10 


Fomalhaut ring (SE ansa) 


5.18 


5.94 


20.26 


3.39 


e Eridani asteroid belt 2 


2.25 


2.07 


6.98 


2.72 


e Eridani sub-mm ring inner edge 


1.52 


1.20 


3.96 


1.69 


e Eridani sub-mm ring outer edge 


0.04 


0.02 


0.11 


0.04 
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Table 6. Fomalhaut 3.6 /xm Dropout Sources 



RA[2000] 


Dec[2000] 


D [AU] 


[4.5] 


[5.8] [8.0] 


Notes 


o^^.oy i ( uu 


-zy.u ( zooi 


ioy i 


1 Q 1 7_i_ n 1 7 
ly . i ( zc u. i ( 






OA A A A -I rAn 


-zy.DDUDUo 


1 n74 

lU / 41: 


1 / .OOlC U.U41: 






OA A AA-\ K/IO 


-zy.DDUDUo 


1 074 

lU / 41: 


1 7 QQ-L n 0/1 

1 / .OOlC U.U41: 






OA A Q70QQQ 


-zy.D41:00 / o 


DoO 


1 7 OR-I- 0/1 

1 / .yoic U.U41: 








90 fiAfi9fi7 

-Zy.D^DZD / 


DDZ 


1 Q 97-1- Ofi 

lo.z / in u.UD 






041:41:. oyyo^z 


-zy.o^i ( 00 


0^0 


in KA4- 00 
iy.04l:ZE U.UU 






O^t^t.^tZ^tOZi 


-zy. 041:0100 


OOZ 


io.OZzC U.iO 






041:41:. OOZoo ( 


9Q fi'ifin/i/i 

-Zy.OOOU4l:4l: 


/191 

41:Z i 


1 o /I Q_i_ n 00 

iO.^OzC U.UU 






344.424175 


-29.644953 


629 


17.29± 0.05 






344.383133 


-29.629925 


278 


17.59± 0.00 






344.419842 


-29.635978 


379 


17.18± 0.07 


17.93± 0.45 




344.427429 


-29.619778 


113 


16.05± 0.00 


14.70± 0.04 14.81± 0.02 


(1) 


344.479287 


-29.629350 


449 


17.59± 0.02 






344.433550 


-29.617836 


177 


18.98± 0.00 






344.400317 


-29.607925 


407 


17.54± 0.07 






344.435804 


-29.613942 


272 


17.80± 0.00 






344.416708 


-29.609828 


348 


18.87± 0.00 






344.421133 


-29.604700 


492 


17.12± 0.04 


15.69± 0.04 


(1) 


344.419792 


-29.600181 


616 


17.89± 0.05 






344.418371 


-29.601414 


581 


18.52± 0.08 






344.413029 


-29.574408 


1329 


18.22± 0.05 


17.25± 0.21 


(1) 


344.423208 


-29.568053 


1506 


17.86± 0.04 


16.22± 0.06 


(1) 



(-'^)[4.5] — [5.8] or [4.5] — [8.0] color suggest source being a background giant 
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Table 7. e Eridani 3.6 //m Dropout Sources 



r> A [9(1(1(1] 

-L V.-i 1- —V./ \.i\.J 


Tlpp [9(10(11 








[■J.OJ 


[8 ni 




53.200488 


-9.512958 


724 


18.32± 


0.02 


16.13± 0.01 


18.20± 0.13 


(1) 


53.194846 


-9.487653 


541 


18.71± 


0.04 




16.45± 0.06 


(1) 


53.210854 


-9.514525 


690 


17.93± 


0.01 


16.79± 0.01 


16.05± 0.00 


(1) 


53.206725 


-9.482892 


402 


18.55± 


0.03 








53.213421 


-9.489692 


417 


18.79± 


0.03 








53.218321 


-9.494831 


448 


18.89± 


0.04 








53.219821 


-9.506739 


574 


18.00± 


0.02 








53.219896 


-9.509417 


604 


18.73± 


0.04 


15.57± 0.01 


16.80± 0.04 


(1) 


53.193787 


-9.440481 


479 


18.49± 


0.03 








53.202425 


-9.455889 


335 


18.64± 


0.03 








53.219771 


-9.481456 


299 


18.38± 


0.04 






(2) 


53.223296 


-9.488964 


366 


18.92± 


0.05 


17.48± 0.06 




(1,2) 


53.213246 


-9.447133 


245 


18.56± 


0.04 






(2) 


53.219871 


-9.466025 


160 


16.35± 


0.01 








53.217404 


-9.469114 


204 


17.21± 


0.02 




15.67± 0.03 


(1) 


53.233725 


-9.477178 


220 


17.55± 


0.02 








53.202654 


-9.420789 


543 


18.10± 


0.02 


. . t 


. . . 




53.209637 


-9.429756 


413 


18.86± 


0.04 








53.241296 


-9.471792 


193 


18.39± 


0.08 








53.239879 


-9.477864 


246 


18.59± 


0.06 








53.209096 


-9.416711 


543 


18.08± 


0.02 








53.225883 


-9.443644 


179 


18.63± 


0.08 








53.243187 


-9.458864 


135 


16.79± 


0.03 


16.57± 0.14 


15.28± 0.02 


(1,3,4) 


53.238017 


-9.463064 


93 


17.11± 


0.10 


16.25± 0.14 


16.42± 0.12 


(1.4) 


53.247821 


-9.466464 


211 


18.07± 


0.04 


17.99± 0.12 


17.09± 0.07 


(1) 


53.258658 


-9.484592 


436 


18.33± 


0.02 






(2) 


53.250783 


-9.458222 


222 


17.31± 


0.02 


16.73± 0.04 


15.69± 0.02 


(1) 


53.229712 


-9.404978 


613 


18.64± 


0.01 








53.245933 


-9.422414 


444 


18.26± 


0.02 








53.256567 


-9.436811 


379 


17.77± 


0.02 








53.255775 


-9.436803 


373 


17.62± 


0.02 
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Table 7 — Continued 



RA[2000] 


Dec[2000] 


D [AU] 


[4.5] 


[5.8] 


[8.0] 


Notes 


53.283862 


-9.494122 


731 


18.97± 0.02 


15.20± 0.00 




(1) 


53.273888 


-9.428683 


595 


17.83± 0.01 


16.35± 0.01 


15.87± 0.00 


(1) 



(1) 



[4.5] — [5.8] or [4.5] — [8.0] color suggest source being a background giant 



Source listed in Table 3. iMarengo et al.l ( 20061 ) 



Source 11 in Table 4, Marengo et all (|2006|) 



Within field of view, but not detected by [Macintosh et al.l (j2003r ) 



A 



Fig. 1. — Correction of the muxstripe artifacts for the first dither position of the 2007 
February 16 e Eridani observations: (a) coadded stack of 256 x 64 subarray images; (h) the 
same stack after PSF subtraction; (c) muxstripe correction matrix; (d) final PSF-subtracted 
stack with correction apphed. The white area at the center of the (b) and (d) panels has 
been masked to exclude pixels where the muxstripe artifacts cannot be efficiently corrected. 
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Fig. 2. — Final roll-subtracted subarray images: (a) Fomalhaut 3.6 fim, (b) Fomalhaut 
4.5 /im, (c) e Eridani 3.6 fim and (d) e Eridani 4.5 fim. The images are scaled linearly from 
—20 to +20 MJy/sr. The ellipse and cross in panels (a) and (b) indicate the position of 
the debris ring and of Fomalhaut b. The dashed circles in panels (c) and (d) indicate the 
inner and outer radii of the e Eridani sub- millimeter ring while the dotted circle shows the 
location of the 20 AU "asteroid belt" (from lBackman et al.ll2009h . The central region of each 
image is masked (because of saturation and PSF subtraction artifacts). The bright spot NE 
of e Eridani is most likely a PSF subtraction artifact, being too narrow (less than 1/2 of 
the PSF FWHM) and lacking the two negative aliases expected for real astronomical point 
sources in the PSF roll-subtracted images. 
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Fig. 3. — Example of PSF roll-subtraction with IRAC full array data. The top panel shows 
the Fomalhaut first epoch image (4.5 /xm). The bottom panel shows the final two-epochs 
roll-subtracted image of Fomalhaut at the same wavelength. Both images are shown in loga- 
rithmic scale, from to 100 MJy/sr in the surface brightness color scale. The cross indicates 
the position of Fomalhaut b. An enlarged figure of the central area around Fomalhaut and 
e Eridani is shown in Section [3^ 
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(a) Fomalhaut 3.6^im 




100 150 
r [arcsec] 



(b) Fomalhaut 4.5/im 



50 100 150 

r [arcsec] 



5 M, 



1 M, 



(c) Fomalhaut 5.8/^m 



5 M, 



^ ^ ~L 




100 150 
r [arcsec] 



(d) Fomalhaut 8.0/im 



V 



. 2 

1 M, 



100 150 
r [arcsec] 



Fig. 4. — 3cr sensitivity of full frame (dashed lines) and subarray (solid lines) Fomalhaut 
images: (a) 3.6 fim, (b) 4.5 fim, (c) 5.8 /im, and (d) 8.0 ^m. The vertical dotted line indicates 
the projected separation of Fomalhaut b. The horizontal lines in dicate the magnitude s 
of 200 Myr p l anets interpolated from 100 and 300 Myr models by iBurrows et al.l (120031 ). 



Fortney et al.l (120081 ) models predict a similar 3 Mj limit at 4.5 fim. 
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100 
T [arcsec] 



(b) e Eridani 4.5/2m 



5 M, 



50 



100 150 
r [arcsec] 



200 



(c) e Eridani 5.8/^m 



10 M, 



y ^ ^ 



V g M, 



50 



100 150 
r [arcsec] 



200 



i? in 



(d) e Eridani B.O^tm 



10 m; 

5 M, 



^ . ~» . 



50 



100 150 
r [arcsec] 



200 



Fig. 5. — 3(7 sensitivity of full frame (dashed lines) and subarray (solid lines) e Eridani 
images: (a) 3.6 /xm, (b) 4.5 fim, (c) 5.8 fim, and (d) 8.0 /xm. The vertical dotted lines 
indicate the inner and outer radii of the submillimeter riri g, and the horizontal lines the 
magnitudes of 1 Gyr planet models by iBurrows et al.l (120031 ). 
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Fig. 6. — Color-magnitude diagram of all point sources detected in the PSF roll-subtracted 
4.5 /im image of Fomalhaut (top) and e Eridani (bottom). Points with error bars are sources 
detected in both 3.6 and 4.5 /im IRAC bands, while the arrows indicate sources that are 
undetected at 3.6 yum (for which we used the image sensitivity at the sou rce location as 



3.6 Atm limit). The squares indicate the colors and magnitudes of planets from lBurrows et al. 



( 120031 ) models, interpolated to the age of each star. The dashed line indicates the average 
3cr sensitivity limit of our 3.6 and 4.5 /xm images. 
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Fig. 7. — Final 4.5 jjm PSF roll-subtracted full frame images. Top panel shows the inner 
50x50" (~ 385x385 AU) of Fomalhaut. The ellipse indicates the location of the debris 
ring, and the cross the position of Fomalhaut b. The bottom panel shows the same area (in 
arcsec, corresponding to ~ 160x160 AU) around e Eridani. The dotted circle indicates the 
position of the 20 AU asteroid belt, while the two da shed circles indi c ate th e inner and outer 
boundaries of the sub-millimeter ring as derived by iBackman et al.l (120091 ) . The images are 
scaled linearly from to 0.3 MJy/sr in the surface brightness color scale. Circle points mark 
the position of 3.6 /im "dropout" sources that are likely background mass-losing giants while 
square points indicate sources detected only at 4.5 /im. 
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Fig. 8. — Foma^ 
mod el used by 



haut b spec t ral en ergy distribution, compared with the iFortney et al.l (l2008l ) 



from Kalas et al. 



Kalas et al.l (120081) to fit their data. Optical photometry and 3cr limits are 



J2008h . Limits are indicated by the tip of the arrows. IRAC 3a limits (thick 
arrows) are from Table HI in this work. The horizontal bars mark the equivalent broad-band 
flux found by integratin g the model spectr um over the instrumental pass-bands. Our 4.5 fim 
limit, together with the lKalas et al.l (120081 ) near-IR limits, suggest a Fomalhaut b mass lower 
than 3 Mj (for a 200 Myr age model), and is inconsistent with the measured optical flux of 
the planet. 



